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DEVELOPING  A  MODEL  FOR  PREDICTING 
SNOWPACK  PARAMETERS  AFFECTING 
VEHICLE  MOBILITY 


Roger  H.  Berger 

INTRODUCTION 

The  primary  mission  of  any  vehicle  is  to  transport  a  payload  from  one  location  to  another.  To 
assess  a  vehicle’s  performance,  we  must  determine  the  constraints  imposed  on  the  vehicle’s  mobility 
by  the  environment  in  which  it  operates.  Because  large  areas  in  the  northern  hemisphere  receive 
seasonal  snowcover  and  because  the  snow  often  severely  limits  vehicle  mobility,  it  is  important  to 
assess  this  problem.  To  do  this,  we  must  determine  how  specific  snow  conditions  affect  vehicle 
performance  and  the  geographical  and  temporal  distribution  of  these  conditions.  This  report 
addresses  part  of  the  latter  problem  by  describing  a  model  that  predicts  the  condition  of  the  snow- 
pack  and  the  snowpack  parameters  that  affect  vehicle  mobility.  The  feasibility  of  this  model  is 
then  examined  in  light  of  the  measurements  that  will  probably  be  available. 

For  studying  mobility  the  snowpack  characteristics  may  be  divided  into  extensive  and  intensive 
properties.  The  more  important  extensive  snowpack  properties  are  depth  and  spatial  and  temporal 
distribution;  the  intensive  properties  are  density,  grain  size  distribution,  temperature  and  free  water 
content.  A  complete  snowpack  model  should  predict  both  classes  of  properties.  To  develop  a 
model  for  predicting  these  quantities,  we  must  understand  the  environmental  factors  that  influence 
the  formation  and  dissipation  of  the  snowpack.  In  general,  the  environmental  effects  with  the 
most  influence  on  the  accumulation  and  the  metamorphosis  and  ablation  processes  in  the  snow¬ 
pack  are  different,  so  these  two  phases  can  be  treated  independently. 

The  accumulation  of  the  snowpack  depends  initially  on  the  snowfall  distribution  within  a  snow¬ 
storm  and  then  on  the  redistribution  by  wind.  Snowstorms  are  of  three  types:  cyclonic,  orographic 
and  convective;  each  has  its  own  pattern  of  snowfall. 

Cyclonic  storms  usually  extend  over  10*-10'^  km^.  The  depth  of  the  snowfall  depends  on  the 
amount  of  moisture  available,  the  strength  of  the  low  and  the  distance  from  the  storm  center.  The 
distribution  of  the  snowfall  is  governed  mainly  by  the  location  of  the  storm  track  and  the  topog¬ 
raphy. 

The  distribution  of  orographic  snowfall  is  governed  by  the  wind  speed  and  direction,  the  eleva¬ 
tion,  the  slope  orientation  and  the  extent  of  shadowing  by  upstream  barriers.  Snowfall  correlates 
well  with  altitude  (Curry  and  Mann  1965,  McPartland  et  al.  1971,  Hendrick  and  DeAngelis  1976). 
The  snow  accumulation  patterns  on  a  scale  of  10^-10*  km^  can  be  predicted  fairly  well  if  the 
direction  of  the  winds  governing  storm  motion  are  known.  The  accumulation  is  sensitive  to  this 
wind  direction  and  can  vary  greatly  with  a  change  in  wind  direction. 

Convective  snowstorms  of  1-10^  km^  can  produce  marked  variations  in  snowfall.  These 
storms  are  generally  responsible  for  the  snowbelts  associated  with  bodies  of  water;  cold  air  passing 


over  warmer  water  produces  snowfall  on  the  lee  side.  This  effect  decreases  with  decreasing  lake 
size  because  of  the  shorter  water  uptake  path  across  small  lakes  and  because  smaller  lakes  freeze 
earlier. 

The  effect  of  large  topographic  barriers  is  well  documented,  using  snowpack  records  from  the 
windward  and  leeward  sides  of  mountain  ranges  (Fiebiger  1976).  The  macroscale  snowpack  defi¬ 
cit  in  the  lee  of  large  barriers  does  not  hold  true  on  a  microscale,  where  the  opposite  frequently 
occurs.  On  this  scale,  wind  redeposits  snow  from  the  windward  side  of  the  crest  of  a  topograph¬ 
ical  feature  to  the  lee  side,  where  it  forms  a  drift  or  cornice.  Vegetation  tends  to  even  out  the 
deposition  of  snow  over  a  crest  by  reducing  or  preventing  wind  transport  of  snow.  The  density, 
type  and  maturity  of  the  vegetation  determines  to  what  extent  the  wind  field  is  modified  (Reif- 
snyder  1955).  The  vegetative  canopy  decreases  the  ground  snowcover  by  intercepting  falling  snow 
(Meiman  1968).  Some  of  this  intercepted  snow  evaporates  or  is  redistributed  under  the  vegetation 
and  in  adjacent  clearings  (FfoUiott  et  al.  1965,  Gary  1975).  As  the  vegetative  cover  increases,  the 
effect  of  the  wind  decreases  until  it  plays  only  a  small  role  in  the  snow  distribution. 

In  alpine  areas  above  the  forest  edge  and  in  open  areas  the  distribution  of  snow  is  controlled  by 
the  wind  and  the  microrelief  of  the  terrain.  Fiebiger  (1976)  found  that  the  reduction  of  wind  in 
forests  produces  a  snow  distribution  nearly  independent  of  the  microrelief.  The  snow  in  forests  is 
generally  deeper  and  less  dense  than  on  open  land,  with  the  difference  correlated  with  the  site,  the 
density  and  the  type  of  trees. 

On  prairies  and  large  open  expanses  the  wind,  the  topography  and  the  vegetation  are  the  domi¬ 
nant  factors  controlling  snow  distribution.  Snow  accumulates  in  hollows  and  streambeds  and  in  the 
lee  of  hills,  ridges  and  clumps  of  trees  or  bushes.  Grasses  and  crop  stubble  can  also  trap  snow  to 
the  height  of  the  vegetative  cover  (Willis  1979). 

The  snowpack  metamorphosis  and  ablation  processes  depend  on  the  net  heat  transfer  between 
the  pack  and  the  environment.  The  processes  may  be  divided  into  two  classes:  radiative  and  non- 
radiative  heat  transfer.  Radiative  heat  transfer  depends  on  the  snow  surface  temperature  and  albedo, 
vegetative  canopy,  cloud  cover  and  insolation  (Anderson  1973,  1976,  Leaf  and  Brink  1973).  Non- 
radiative  heat  transfer  is  through  conduction,  convection,  latent  heat  exchange  and  mass  exchange. 
These  processes  are  functions  of  the  snow,  air  and  ground  temperatures,  the  dew  point  and  the 
wind  velocity  (Leaf  and  Brink  1973,  Anderson  1976). 


REVIEW  OF  EXISTING  MODELS 

Most  of  the  extensive  body  of  literature  on  snowcover  modeling  has  been  developed  for  hydro- 
logical  purposes.  For  the  hydrologist  the  melting  or  ablation  phase  of  the  snowpack  is  of  primary 
interest,  and  accumulation  processes  are  only  of  secondary  interest.  The  general  hydrological 
approach  is  to  develop  a  model  of  the  snowpack  based  on  energy  considerations.  This  model  is 
calibrated  against  the  meteorological  and  hydrological  records  of  a  specific  area  or  watershed  to 
determine  the  effects  of  various  unmeasured  parameters.  The  model  is  then  used  to  predict  runoff. 

The  basic  understanding  of  the  relationships  between  the  environment  and  the  behavior  of  the 
snowpack  that  are  used  in  most  of  the  current  models  can  be  traced  to  the  pioneering  work  de¬ 
scribed  in  Snow  Hydrology  (U.S.  Army  Corps  of  Engineers  1956).  This  publication  describes 
studies  of  three  western  mountain  watersheds  and  the  application  of  the  basic  sciences  to  the 
hydrological  cycle.  In  Snow  Hydrology  each  component  of  the  hydrological  cycle  was  determined 
for  the  given  environmental  conditions.  To  do  this,  the  mathematical  relationships  between  the 
variables  affecting  each  component  were  derived  and  constants  of  proportionality  were  determ¬ 
ined.  With  these  relationships  the  behavior  of  the  snowpack  at  a  point  was  determined;  these 
point  values  were  then  used  to  determine  the  distribution  and  behavior  of  the  snowpack  with  re¬ 
spect  to  environmental  differences. 
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For  hydrological  purposes  it  is  impractical  and  unnecessary  to  use  the  kind  of  point-by-point 
analysis  used  in  Snow  Hydrology.  A  more  practical  approach  is  1)  to  determine  averages  for  each 
major  topographical  and  environmental  subdivision  within  a  basin  and  2)  to  ignore  individual  events 
and  deal  with  time  averages.  This  concept  leads  to  the  use  of  indexes  to  represent  basin  averages. 

To  use  indexes  properly,  the  physical  character  of  the  area  and  of  each  measurement  point  must  be 
considered.  Since  the  use  of  indexes  involves  sampling  theory  and  measurement  errors,  statistical 
procedures  may  be  used  to  establish  the  reliability  of  estimates  and  the  weightings  of  individual 
factors.  Indirect  methods  can  then  be  used  to  evaluate  basin  parameters,  using  observational  data 
commonly  available;  however,  these  methods  must  be  used  with  care,  and  additional  data  must  be 
gathered  when  necessary. 

The  conclusions  reached  in  Snow  Hydrology  were  corroborated  by  a  similar  Russian  study. 
Melting  of  Snow  Cover  (Kuz'min  1961).  While  the  two  investigations  overlap,  the  Russian  study 
also  includes  methods  that  apply  to  plains  regions,  where  the  snowcover  tends  to  be  shallow  and 
discontinuous  rather  than  deep  and  continuous,  as  in  most  mountain  regions. 

The  snowpack  modeling  that  is  the  direct  descendant  of  Snow  Hydrology  and  Melting  of  Snow 
Cover  can  be  divided  into  three  types;  1)  energy  and  mass  balance  models,  2)  radiation  balance  and 
energy  exchange  models,  and  3)  temperature  index  models.  These  model  types  are  all  based  on  the 
same  heat  transfer  method  or  description  of  the  snowpack;  the  difference  is  in  the  detail  with  which 
the  physical  processes  within  the  snowpack  and  the  interactions  of  the  environmental  factors  and 
the  snowpack  are  treated.  None  of  these  models  treat  the  accumulation  process  in  detail. 

The  models  that  treat  the  accumulation  process  are  more  diverse  in  their  underlying  principles. 
They  are  1)  orographic  precipitation  models  and  2)  statistical  snowfall  models.  The  statistical 
models  are  based  on  the  statistical  treatment  of  either  the  meteorological  processes  governing  snow¬ 
fall  or  historical  snowfall  records. 

Accumulation  models 

Orographic  precipitation  models 

These  models  are  based  on  the  observation  that  precipitation  increases  with  increasing  elevation 
in  mountainous  regions.  However,  linear  regressions  between  elevation  and  precipitation  (Storr  and 
Golding  1974,  Curry  and  Mann  1965)  yield  relationships  that  apply  only  to  the  localities  for  which 
they  were  developed.  This  suggests  that  elevation  is  not  the  only  terrain  variable  affecting  precip¬ 
itation  in  mountainous  regions.  Spreen  (1947)  found  that  30%  of  the  variation  in  precipitation 
could  be  accounted  for  by  elevation,  but  by  including  terrain  slope,  exposure  and  orientation  in 
addition  to  elevation,  88%  of  the  variation  could  be  explained.  Terrain  slope  is  the  maximum 
change  in  elevation  within  a  five-mile  radius,  exposure  is  a  measure  of  the  number  of  barriers  with¬ 
in  a  20-niile  radius  that  have  elevations  exceeding  that  of  the  point  of  interest  by  more  than  1000 
ft,  and  orientation  is  the  direction  of  greatest  exposure. 

In  mountainous  terrain  the  orographic  precipitation  process  is  dominant.  This  process  can  be 
expressed  by  the  following  formula  (Rhea  1978a): 

«=  i  f 

0»*  ./O  Jp, 

where 

R  =  precipitation 
E  =  precipitation  efficiency 
p^  =  density  of  water 
g  =  acceleration  due  to  gravity 


dz 


wdPdt 


(1) 
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f. 


dt  =  time  interval 


rP2 

Jp  dP  =  vertical  extent  of  the  saturated  air  mass  producing  the  precipitation  (expressed  in 
*  terms  of  pressure). 


=  vertical  gradient  of  specific  humidity 
w  =  vertical  velocity  of  the  air  parcel. 


In  orographic  precipitation  the  vertical  motion  in  eq  1  is  produced  by  the  slope  of  the  terrain.  The 
upward  motion  of  the  moist  air  causes  condensation,  some  of  which  reaches  the  ground  as  precipi¬ 
tation. 

Orographic  models  couple  this  precipitation  process  with  slope,  exposure,  elevation,  upper  air 
moisture  and  wind  to  predict  the  distribution  of  snow  accumulation.  Rhea  (1978a)  used  a  two- 
dimensional  air  flow  model,  in  which  the  topography  forces  vertical  displacements  of  the  air  but 
does  not  allow  the  air  to  deviate  horizontally  around  topographic  barriers.  The  model  keeps  track 
of  the  condensation  o:  evaporation  resulting  from  the  vertical  motions  of  the  air  at  the  mesh  points 
of  a  square  topographic  grid.  The  model  is  capable  of  using  a  grid  as  small  as  2.5  km  over  an  area 
of  10*  km^. 

Orographic  models  are  designed  to  use  current  weather  information  to  predict  the  water  resources 
stored  in  the  snowpack.  They  can  yield  realistic  snow  accumulations  averaged  over  a  season,  using 
as  input  the  atmospheric  profiles  from  a  few  widely  spaced  (300-400  km)  measurement  points. 

The  predictions  are  often  inaccurate  for  particular  storms  because  the  distribution  of  precipitation 
is  based  on  the  mesoscale  wind  direction  and  doesn’t  allow  for  the  horizontal  wind  direction  changes 
due  to  local  terrain  features.  These  models  also  ignore  cyclonic  storms,  which  are  not  significant  on 
a  seasonal  basis  but  can  be  very  significant  on  a  short-term  basis. 

For  mobility  studies,  predictions  for  individual  storms  must  be  more  accurate.  This  may  be 
accomplished  by  scaling  down  the  size  of  the  region  of  prediction  and  by  using  a  three-dimensional 
flow  model  with  more  detailed  wind  input. 

Statistical  snowfall  models 

The  first  of  these  models  was  developed  by  Sage  (1976).  Its  object  is  to  produce  a  synthetic 
hourly  snowfall  sequence  that  is  statistically  indistinguishable  from  a  real  snowfall.  The  model  uses 
partitioning,  in  which  each  hour  of  a  storm  is  designated  as  a  period  with  or  without  snowfall.  The 
sequence  of  snowfall  periods  is  then  modeled  using  a  first-order  regression  equation.  Similar  equa¬ 
tions  are  developed  for  sequences  of  fall  periods  and  snowfall  intensity  The  three  sequences  are 
then  combined  to  represent  the  hourly  sequence  of  snowfall  intensities. 

The  second  kind  of  model  mathematically  describes  the  statistics  of  individual  snowfall  events 
and  the  meteorological  factors  affecting  the  snowfall.  The  increase  in  the  snowpack  due  to  an  in¬ 
dividual  storm  can  be  expressed  as  (Isyumov  and  Davenport  1974,  Isyumov  and  Mikitiuk  1976); 

AR=gis,p^,T,U)  (2) 


where 

s  =  snowfall  depth 
Pj  =  snow  density 
T=  air  temperature 
U  =  wind  speed. 
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Since  the  meteorological  variables  are  random,  AR  is  also  a  random  variable. 

The  probability  that  a  snowfall  event  will  occur  depends  on  the  probability  that  the  associated 
meteorological  factors  occur.  By  examining  the  interdependence  of  the  meteorological  factors, 
the  principal  dependencies  can  be  identified.  The  probability  that  the  meteorological  conditions 
exist  to  produce  AR  can  be  expressed  as 

p(s,  P5,  T,U)^  PipjT)p(slT)p(slU)p(s)  (3) 

where  p(j)JT)  is  the  probability  density  function  for  p^,  given  the  occurrence  of  T;  this  expresses 
the  fact  that  Tis  the  primary  meteorological  variable  on  which  Pj  depends.  These  joint  probabil¬ 
ities  are  determined  empirically. 

A  snowfall  of  depth  s  is  a  random  event.  The  snowfall  process  is  not  stationary  but  varies  with 
time  during  the  snow  season.  Isyumov  and  Davenport  (1974)  showed  that  the  snowfall  statistics 
on  a  24-hour  basis  fit  a  Weibull  distribution,  and  thus  the  snowfall  probability  may  be  expressed  as 

p(>s,  t)  =  m(t)  exp|-  |  (4) 

where 

p(>j,  /)  =  conditional  probability  that  a  snowfall  of  depth  greater  than  s  occurred  at  time  t 
m(t)  =  snowfall  expectation  at  time  t 
C(t),K(t)  =  WeibuU  distribution  parameters. 

The  three  time  functions  are  obtained  by  fitting  the  snowfall  distribution  for  the  entire  season  to  a 
Weibull  distribution.  These  parameters  are  measures  of  the  length  of  the  snow  season  and  the  rela¬ 
tive  distribution  of  total  snowfall  nt(t),  the  probability  of  large  snowfalls  K(t),  and  the  mode  of  the 
snowfall  magnitude  distribution  C(/). 

Expressions  similar  to  eq  4  are  developed  in  an  analogous  manner  for  the  other  meteorological 
factors  and  are  combined  in  eq  3.  This  equation  is  used  in  a  Monte  Carlo  calculation  to  simulate 
the  day-by-day  snowfall. 

For  mobility  studies  the  short-range  predictions  of  statistical  snowfall  models  are  not  adequate, 
since  they  have  no  provision  for  including  current  meteorological  data.  The  utility  of  these  models 
may  lie  in  Monte  Carlo  simulations  for  mobility  modeling,  rather  than  in  actual  snowpack  predic¬ 
tions. 

Ablation  models 

Since  the  temperature  index,  the  radiation  balance  and  energy  exchange,  and  the  energy  and 
mass  balance  snowpack  models  are  all  based  on  the  same  physical  description  of  the  snowpack 
energy  exchange  processes,  the  most  detailed  will  be  described  first,  followed  by  a  description  of 
the  simplifications  used  in  the  other  models. 

Energy  and  mass  balance  models 

These  models  describe  the  ablation  and  metamorphosis  of  the  snowpack  in  terms  of  two  equations, 
the  energy  balance  equation  and  the  equation  for  heat  transfer  within  the  snowpack.  The  snowpack 
is  divided  into  layers,  and  implicit  finite  difference  equations  are  written  for  each  layer.  This  set 
of  equations  is  then  solved  by  an  iteration  technique.  Included  in  the  models  are  the  densification 
of  the  snow,  the  change  in  thickness  of  each  layer,  and  the  retention  and  transmission  of  water. 

The  energy  balance  equation  is  derived  from  the  snowpack  energy  inputs  and  losses  that  occur 
at  the  air/snow  and  soil/snow  interfaces.  Under  most  conditions  the  energy  exchange  at  the  air/ 
snow  interface  is  much  greater  than  at  the  soil/snow  interface.  The  energy  exchanges  across  the 
interfaces  are  either  by  the  transfer  of  mass  or  by  radiation.  The  mass  transfer  is  in  the  form  of 


water  in  any  of  its  three  phases.  Water  is  lost  by  evaporation  and  sublimation  at  the  air/snow  in¬ 
terface  and  liquid  water  outflow  at  the  soil/snow  interface.  Water  flows  from  the  bottom  of  the 
snow  cover  when  the  snow  is  isothermal  at  0°C  and  saturated  with  water,  when  the  snow  is  referred 
to  as  “ripe.”  The  quantity  of  interest  in  these  instances  is  the  transfer  of  latent  or  sensible  heat. 

The  net  energy  transfer  by  radiation  is  the  difference  between  the  incoming  and  outgoing  radia¬ 
tion.  Both  visible  and  infrared  wavelengths  have  to  be  considered  because  a  significant  amount  of 
energy  can  be  absorbed  in  the  visible  and  reradiated  in  the  infrared.  The  incoming  radiation  comes 
directly  from  the  sun  or  is  reflected  and/or  emitted  by  the  atmosphere  and  terrestrial  objects.  The 
outgoing  radiation  is  made  up  of  reflected  solar,  reflected  long-wave,  and  emitted  radiation.  For 
most  modeling  purposes  the  emissivity  e  of  snow  in  the  infrared  may  be  taken  as  one,  and  the  in¬ 
frared  reflectivity  and  emissivity  are  related  by  e  =  1  -  p.  The  reflected  long-wave  component  may 
be  neglected.  The  net  radiative  transfer  can  be  expressed  as  (Anderson  1976) 

Q„  =  (2i-(2,  +  (6-(2a)-(A/-6-a-7:^)  (5) 

where 

Qj  =  insolation 

Qj  =  reflected  solar  radiation 
e-Qj  =  incident  atmospheric  and  terrestrial  radiation 
A/-e-a-7’^=  radiation  emitted  by  the  snowpack  with  surface  temperature  Tq  during  computation 
interval  At 

a  =  Stefan-Boltzman  constant. 

Several  environmental  factors  affect  the  net  radiative  energy  transfer.  The  incident  radiation 
depenas  on  factors  such  as  latitude,  atmospheric  clarity  and  cloud  cover,  topography  and  vegeta¬ 
tion  cover.  The  outgoing  radiation  is  a  function  of  the  albedo  and  the  surface  temperature  of  the 
snowcover;  both  vary  on  a  small  scale,  but  the  surface  temperature  is  also  influenced  by  larger 
scale  variations. 

The  nonradiative  heat  transfer  may  be  considered  as  four  processes,  the  first  three  occurring  at 
the  air/snow  boundary:  1)  sensible  heat  transfer,  2)  latent  heat  transfer,  3)  transfer  by  mass 
changes,  and  4)  one  general  heat  transfer  process  at  the  soil/snow  boundary.  The  sensible  and 
latent  heat  transfer  depends  on  the  condition  of  the  atmospheric  boundary  layer,  which  is  turbu¬ 
lent  under  most  conditions.  Within  this  boundary  layer  there  is  a  thin  {~  1  mm)  viscous  sublayer 
where  molecular  conduction  processes  are  important.  Above  this  viscous  layer  is  a  layer  (meters  to 
tens  of  meters  thick)  where  shear  stress  and  vertical  fluxes  of  heat  and  water  vapor  can  be  considered 
independent  of  height.  It  is  this  latter  layer  that  is  used  in  estimating  latent  and  sensible  heat  exchange. 

If  the  hypothesis  is  correct  that  the  mean  vertical  fluxes  of  latent  and  sensible  heat  are  propor¬ 
tional  to  their  vertical  gradients,  then  the  water  vapor  transfer  is  predominantly  a  function  of  the 
wind  speed  and  the  vapor  pressure  gradient.  The  particular  relationship  depends  upon  the  lapse 
rate  and  turbulence,  as  characterized  by  the  Richardson  number. 

Since  a  full  set  of  measurements  at  two  heights  above  the  snow  surface  is  not  practical  in  most 
cases,  the  expressions  for  the  latent  heat  transfer  and  the  sensible  heat  transfer  are  given  in 
terms  of  the  snow  surface  quantities  and  turbulent  layer  measurements: 

(6) 

Gh=^s- ^  •‘y/W-(7’.-7’o)  (7) 
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where 


=  latent  heat  of  sublimation 
=  density  of  water 
f{U)  =  wind  speed  function 

^0  “  vapor  pressure  in  the  air  and  at  the  snow  surface,  respectively 
7  =  psychrometric  constant 

T^,  Tq  =  temperature  of  the  air  and  the  snow  surface,  respectively. 

The  snowpack  also  gains  and  loses  heat  by  mass  changes  resulting  from  precipitation  and  vapor 
transfer.  If  the  heat  content  of  the  vapor  is  assumed  to  be  negligible  and  the  water  outflow  is 
assumed  to  be  at  0°C,  only  the  heat  transferred  by  precipitation  must  be  considered.  The  precipi¬ 
tation  temperature  is  assumed  to  be  the  wet-bulb  temperature,  and  the  heat  transfer  is  expressed 
as 

Gm=C^Pw'R(7’w-273.16)/10  (8) 

where 

=  specific  heat  of  water 
p^  =  density  of  water 
R  =  precipitation  (mm) 

=  wet-bulb  temperature  (/T)  . 

These  four  terms  represent  the  energy  exchange  at  the  snow/air  interface.  An  additional  term 
is  needed  for  the  exchange  at  the  snow/soil  interface.  The  sum  of  these  five  terms  is  equal  to  the 
change  in  snowpack  energy  content: 

A0  =  (x.p,)'  (9) 

where 

X  =  snow  depth 
C.  =  specific  heat  of  ice 

T,  =  temperature  of  the  ice  grains  inside  the  snowpack 
Lf  =  latent  heat  of  fusion 

W  =  amount  of  liquid  water  in  the  snow,  expressed  as  depth 
^  =  time  since  the  last  computation. 

The  complete  energy  balance  equation  is; 

(2„+<2e+Gh+Gm+Gg-^G  =  o  (10) 

where  is  the  heat  transfer  at  the  surface  of  the  ground.  Combining  eqs  5-9  gives 

[G,-Gr  +  (e-G.)-(Af-e-a.7^^  +  [l,  •  ^  /(G)  •  (e.  -  eo)| 

+  •  7  •  (7-.  -  T-o)]  +  [c„Pw«  {T^  -  273.16)/lo] 

+  Gg- j(^  P.)'  [(C,-7’.)'"^'-(C,T/]  -  ^  (K/'-Af_^)|=  0  (11) 

This  is  the  first  equation  on  which  the  model  is  based;  the  second  expresses  the  heat  transfer 
within  the  snowpack.  Heat  is  transferred  within  the  snowpack  by  conduction,  radiation,  phase 


change  and  vapor  diffusion.  The  mathematical  expressions  for  these  mechanisms  include  coef¬ 
ficients  that  are  related  to  the  density,  grain  size  and  crystal  structure  of  the  snowpack.  Since 
the  latter  two  quantities  are  not  generally  measured,  the  density  is  the  only  parameter  used  to 
characterize  the  snowpack. 

The  heat  transfer  within  the  snowpack  is  essentially  one-dimensional,  with  heat  either  gained  or 
lost  through  the  air/snow  and  soil/snow  interfaces.  The  one-dimensional  Fourier  heat  conduction 
equation  is  used,  with  modifications  to  include  the  different  mechanisms  of  heat  transfer; 


(12) 

where 

Cj  =  specific  heat  of  ice 
C  =  concentration  of  water  vapor 
=  effective  diffusion  coefficient 
/  =  short-wave  radiation  flux 
Kg  =  effective  thermal  conductivity 
Lj  =  latent  heat  of  sublimation 
Lf  =  latent  heat  of  fusion 
T  =  snow  temperature  . 

The  terms  in  this  equation  account  for  conduction,  heat  sources  and  sinks,  and  diffusion. 

The  one-dimensional  terms  accounting  for  the  heat  conduction  within  the  snowpack  are 
CjPj  3773/  -  3(Kg  37’/3z)/3z.  The  effective  conductivity  is  the  sum  of  the  conduction  by  the  ice  in 
the  snow  fabric  and  by  the  air  in  the  pores,  and  the  radiant  exchange  across  the  pores.  The  second 
term  also  accounts  for  the  fact  that  this  effective  conductivity  is  a  function  of  both  the  position 
in  the  snowpack  due  to  layering  and  the  change  in  density  with  depth.  Heat  sources  and  sinks  are 
accounted  for  by  Lf  3H^/3f  -  3//3z.  Water  within  the  snowpack  acts  as  either  a  source  or  a  sink  of 
heat,  depending  upon  whether  it  is  decreasing  due  to  freezing  or  increasing  due  to  melting.  These 
effects  are  accounted  for  by  the  first  term;  the  second  term  represents  the  heat  due  to  insolation. 
Since  the  short-wave  solar  radiation  penetrates  to  significant  depths  before  it  is  absorbed,  it  acts 
as  an  extended  heat  source  in  the  upper  layers  of  the  pack. 

Temperature  gradients  within  the  snow  are  the  main  driving  forces  behind  the  diffusion  of  water 
vapor  in  the  snowpack.  Associated  with  the  diffusion  process  are  sublimation  and  condensation, 
which  are  the  source  and  sink  for  the  water  vapor,  respectively.  Heat  transfer  due  to  the  transport 
of  the  vapor  is  negligible,  but  the  heat  absorbed  and  released  in  sublimation  and  condensation  is 
significant.  These  processes  decrease  the  internal  temperature  gradient.  These  processes  and  the 
metamorphic  processes  within  the  snowpack  are  accounted  for  in  the  terms  [3C/37’  •  3(Z)g 
37’/3z)/3z  +  Dg  32C/37^  (37’/3z)2] . 

The  effective  conductivity  is  a  function  of  the  distribution,  size,  shape  and  intergranular  bond¬ 
ing  of  the  ice  crystals  but  it  can  be  empirically  correlated  with  density  alone.  The  effective  diffu¬ 
sion  coefficient  can  be  related  to  the  snowpack  temperature,  and  the  absorption  coefficient  can  be 
related  to  density  if  the  local  relationship  between  grain  size  and  density  is  empirically  established. 

The  energy  and  mass  balance  model  is  based  on  the  combined  solution  of  eqs  1 1  and  12. 

The  information  needed  to  solve  the  equations  is  air  temperature,  air  vapor  pressure,  wind  speed, 
incoming  solar  radiation,  reflected  solar  radiation,  incoming  long-wave  radiation,  and  soil  temper¬ 
ature.  This  extensive  set  of  data  restricts  the  use  of  this  model  to  a  few  well-instrumented  water¬ 
sheds. 


Radiation  balance  and  energy  exchange  models 

These  models  use  a  less  complete  description  of  the  energy  exchange  and  heat  transfer  processes 
in  the  snowpack  than  that  used  in  the  energy  and  mass  balance  models.  The  models  compute  the 
energy  exchange  by  one  set  of  equations  and  the  heat  deficit  of  the  snowpack  by  another  equation. 
For  example,  Leaf  and  Brink  (1973)  use  the  air  temperature  and  the  incident  solar  radiation  to 
estimate  the  net  radiation  balance,  from  which  the  energy  exchange  is  computed.  The  one-dimen¬ 
sional  Fourier  heat-conduction  equation  is  used  to  compute  the  temperature  at  the  midpoint  of  the 
snowpack,  assuming  that  the  snow  surface  temperature  equals  the  air  temperature.  The  internal 
and  surface  temperatures  are  used  to  estimate  the  heat  deficit  in  the  snowpack.  This  simplified 
treatment  reduces  the  input  data  requirements  while  producing  a  sufficiently  accurate  snowpack 
simulation. 

Temperature  index  models 

These  models  use  the  air  temperature  as  the  sole  index  of  snowpack  ablation.  Maximum  and 
minimum  temperatures  and  degree-days  are  used  most  commonly.  Anderson  (1973)  developed  a 
heat-transfer  equation  that  is  a  function  of  air  temperature,  so  that  the  heat  exchange  at  the  snow¬ 
pack  surface  can  be  related  to  the  air  alone.  The  simplifications  and  assumptions  involved  lead  to 
a  loss  of  precision,  and  knowledge  of  the  internal  metamorphic  processes  is  lost.  The  big  advantage 
of  the  procedure  is  that  much  less  input  data  is  required  without  losing  much  prediction  accuracy. 

Using  existing  models  for  studying  mobility 

The  existing  models  of  snowpack  accumulation  and  ablation  would  have  to  be  modified  sub¬ 
stantially  before  they  could  be  used  for  mobility  modeling  or  prediction.  Of  the  accumulation 
models,  only  the  orographic  model  appears  to  be  adaptable  for  studying  mobility.  This  could  be 
accomplished  by  reducing  the  grid  size  and  expanding  the  two-dimensional  wind  flow  analysis  to 
three  dimensions.  Unfortunately  this  model  is  restricted  to  areas  of  relatively  large  topographical 
relief,  and  atmospheric  moisture  and  wind  profiles  are  needed.  The  other  two  accumulation  mod¬ 
els  may  be  of  some  utility  in  vehicle  systems  evaluation  modeling  but  are  not  suitable  for  predict¬ 
ing  actual  snow  conditions. 

The  best  ablation  model  for  studying  mobility  will  be  one  incorporating  a  compromise  between 
the  amount  of  input  data  required  and  the  need  for  snowpack  detail.  The  most  reasonable  approach 
is  to  use  a  temperature  index  model  and  to  incorporate  as  much  snowpack  detail  as  possible. 


PROPOSED  SNOWPACK  MODEL  FOR  MOBILITY  STUDIES 

The  existing  snowpack  models  are  capable  of  predicting  the  depth,  density  and  temperature  of 
the  snowpack.  However,  they  are  limited  because  they  cannot  predict  accumulation  events,  they 
can  predict  snowfall  distribution  only  for  Umited  areas,  they  require  extensive  data,  and  they  can- 
no_t  account  for  wind  effects. 

'fhe  models  that  can  predict  many  of  the  snowpack  parameters  needed  for  studying  mobility 
are  only  o*"  limited  use  because  they  require  a  knowledge  of  air  temperature,  dew  point,  wind 
speed,  incoming  and  reflected  solar  radiation,  and  atmospheric  long-wave  radiation.  To  extend 
the  models  from  point  predictors  to  areal  predictors  the  slope,  aspect,  vegetative  cover,  snowpack 
albedo,  and  wind  field  must  also  be  estimated. 

Considering  this  set  of  input  requirements  and  the  errors  involved  in  estimating  the  areal  values 
and  comparing  the  accuracy  of  predictions  of  the  conceptual  models  with  an  index  model,  such 
as  a  temperature  index  model,  the  advantages  of  an  index  model  are  apparent  for  mobility  pur¬ 
poses.  In  addition,  when  the  probable  input  data  set  (which  will  be  available  under  operational 
conditions)  is  considered,  the  index  model  becomes  more  attractive.  The  operational  data  set  is 
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assumed  to  include  the  data  from  three  sources:  1)  a  master  data  bank  of  geographical  and  terrain 
information,  2)  reconnaissance  of  the  snow  and  terrain  conditions,  and  3)  weather  forecasts. 

An  existing  model  that  can  operate  from  this  data  base  has  been  selected  as  the  basis  for  the 
mobility  snowpack  model;  it  is  the  Snow  Accumulation  and  Ablation  Model  (Anderson  1973) 
developed  for  the  National  Weather  Service  River  Forecast  System.  This  model  indexes  each  phys¬ 
ical  process  that  occurs  in  the  snowpack  to  the  air  temperature.  Although  this  results  in  a  signifi¬ 
cant  loss  in  the  accuracy  with  which  the  effects  of  each  process  is  known,  this  is  compensated  for 
by  the  great  reduction  in  input  data  and  the  applicability  to  a  much  greater  geographical  area  where 
only  sparse  data  are  available.  Anderson  (1976)  has  shown  that  his  model  can  predict  about  as 
accurately  as  a  detailed  energy  and  mass  balance  model.  The  agreement  of  the  two  models  is  slightly 
misleading  because  a  more  extensive  set  of  parameters  must  be  determined  for  the  temperature  in¬ 
dex  model  before  it  can  be  used.  The  development  of  a  method  for  determining  these  parameters 
from  weather  station  data  and  extrapolation  techniques  is  one  of  the  areas  where  further  develop¬ 
ment  is  needed. 

Model  components 

In  the  model  the  snowpack  is  considered  to  be  a  dynamic  energy  reservoir,  the  energy  being  in 
the  form  of  heat.  The  capacity  of  the  reservoir  depends  on  the  masses  of  ice  and  water  in  the  snow 
and  their  temperatures.  Changes  in  energy  content  are  measured  relative  to  a  dry  isothermal  snow¬ 
pack  at  0°C.  Heat  lost  by  such  a  snowpack  causes  an  internal  temperature  drop  as  the  energy  is 
conducted  to  the  surface,  where  it  is  lost  to  the  environment.  Heat  added  to  the  pack  will  reduce 
any  energy  deficit  and  produce  mel**ng  in  those  parts  of  the  pack  that  reach  0°C.  These  energy 
exchanges  between  the  snowpack  and  the  environment  are  expressed  in  the  energy  balance  equa¬ 
tion,  which  is  the  first  major  component  of  the  model  and  is  described  below.  The  second  major 
component  is  the  heat  transfer  equation,  which  deals  with  the  temperature  structure  and  heat 
transfer  within  the  snowpack. 

Heat  balance  equation 

The  heat  balance  equation  relates  the  net  change  in  the  energy  content  of  the  snowpack  to  the 
algebraic  sum  of  the  heat  gained  from  or  lost  to  the  environment  by  the  various  mechanisms  of 
heat  transfer: 

AG  =  G„  +  G.  +  G,  +  G„+Gg.  (13) 

All  of  the  terms  except  Gg  represent  processes  of  heat  transfer  across  the  snow/air  interface. 

The  radiative  heat  transfer  is  the  sum  of  the  incident  short-wave  solar  radiation,  the  long-wave 
radiation  from  the  sky  and  forest  cover,  and  the  long-wave  radiation  emitted  by  the  snow.  The 
short-wave  component  of  the  heat  balance  is  affected  by  the  slope,  aspect  and  latitude  of  the  site, 
the  cloud  cover,  the  vegetative  cover,  and  the  short-wave  reflectivity  of  the  snow.  The  long-wave 
components  are  affected  by  the  cloud  cover,  the  vegetative  cover,  and  the  emissivities  and  temper¬ 
atures  of  the  various  emitting  objects.  Gn  expressed  as 

Qn  =  (1-w)  /+  C,  o{l-CJ  ae,  7^  +  C^oc^tJ 

-C^ae^T^  -(l-C^)aeo7^  (14) 

where 

Cj,  Cy  =  cloud  and  vegetative  cover  factors,  respectively 
S,A,L  =  slope,  aspect  and  latitude  factors,  respectively 
CO  =  snow  albedo 
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/=  incident  short-wave  energy 

a  =  empirically  determined  constant  (0.757)  (U.S.  Army  Corps  of  Engineers 
1956) 

Cj,  Cy,  Cq  =  emissivity  of  the  sky,  the  foliage  and  the  snow,  respectively 
T",,  Ty,  Tq  =  temperature  of  the  air,  the  foliage  and  the  snow,  respectively. 

The  first  term  in  eq  14  is  the  net  incident  short-wave  energy.  The  next  two  terms  are  the  incident 
long-wave  radiation  from  the  sky  and  the  foliage,  respectively.  The  remaining  two  terms  are  the 
outgoing  long-wave  energy  to  the  sky  and  the  foliage,  respectively. 

Combining  the  snow-canopy  exchange  terms  and  the  snow-sky  exchange  terms  yields 

=  C,  (I-CJ SA-L  (l-co) I+aiT^C,a-  7^)  (l-CJ  +  (7^  -  7^)  (15) 


where  the  emissivities  are  assumed  to  be  equal  to  one. 

When  the  energy  emitted  by  the  snow  is  calculated,  the  snow  surface  temperature  is  assumed  to 
be  equal  to  the  air  temperature  for  temperatures  less  than  0°C  and  equal  to  0°C  for  air  tempera¬ 
tures  above  0°C.  When  the  ambient  temperature  is  below  0°C  the  net  radiation  exchange  between 
the  canopy  and  the  snow  is  assumed  to  be  zero,  assuming  that  7*0  =  Ty  During  periods  of  precipi¬ 
tation  the  net  long-wave  balance  is  assumed  to  be  zero. 

The  latent  heat  transfer  is  due  to  the  transfer  of  water  vapor  between  the  snow  and  the  air. 
This  process  is  governed  by  the  vapor  pressure  gradient  and  the  wind  speed  and  can  be  expressed 
by  a  relationship  derived  by  Dalton: 

K=/(C0(e,-eo)  0^) 

w^iere 

V  =  water  vapor  transfer 
f(U)  -  wind  speed  function 

c,,  Cq  =  vapor  pressure  of  the  air  and  the  snow  surface,  respectively. 

The  wind  function  is  usually  represented  as 

f(,U)  =  a  +  bU  (17) 

w^iere  U  is  the  wind  speed  and  a  and  b  are  empirically  determined  constants. 

The  values  of  the  constants  depend  on  the  lapse  rate  and  the  stability  of  the  surface  layer.  Since 
stable  conditions  predominate  over  snow,  the  amount  of  vapor  transfer  in  the  absence  of  wind 
(represented  by  the  constant  a)  should  be  small  or  zero  (U.S.  Army  Corps  of  Engineers  1956). 

The  latent  heat  transfer  is  the  amount  of  heat  gained  by  the  snow  through  condensation  lost 
through  sublimation.  This  heat  is  expressed  by 


The  sensible  heat  transfer  Qf,  can  be  obtained  from  the  ratio  of  the  sensible  heat  transfer  to  the 
latent  heat  transfer  known  as  Bowen’s  ratio  (Anderson  1973); 


where  y  is  the  psychrometric  constant. 


On  substitution  of  the  latent  heat  transfer  the  sensible  heat  transfer  becomes 

Qu=bUp^L,y(.T,-T^).  (20) 

The  precipitation  input  to  the  model  will  be  a  forecast  with  both  the  type  and  amount  of  pre¬ 
cipitation  given.  For  the  heat  transfer  calculation  the  temperature  of  rain  is  assumed  to  be  the 
same  as  the  air  temperature.  The  temperature  of  snow  is  taken  to  be  0°C  at  air  temperatures  above 
freezing  and  the  same  as  the  air  temperature  if  it  is  below  freezing.  The  heat  gain  from  rain  is 
applied  against  any  net  heat  deficit  in  the  pack,  and  any  excess  heat  after  the  deficit  is  made  up 
goes  into  melting  the  pack.  In  the  case  of  snow,  the  heat  transfer  produces  no  change  or  increases 
the  net  energy  deficit.  This  heat  transfer  is  expressed  as 

Q„=/?Cp(7’p-7o)  (21) 

where  is  the  specific  heat  of  water  or  ice,  depending  on  the  precipitation  form  (either  C,  or  C^), 
and  Tp  and  Tq  are  the  precipitation  and  snowpack  temperatures,  respectively. 

The  heat  transfer  between  the  snowpack  and  the  ground  is  due  mainly  to  conduction.  Since 
the  ground  can  be  considered  to  be  a  semi-infinite  heat  reservoir,  its  temperature  changes  only  very 
slowly.  Under  these  conditions  a  steady-state  heat  conduction  equation  can  be  used  to  express  the 
heat  transfer: 

(22) 

where  Kg  is  the  conductivity  of  the  soil  and  bT/dz  is  the  temperature  gradient  in  the  soil  at  the 
soil/snow  interface. 

The  total  heat  transfer  is  obtained  by  substituting eqs  IS,  18, 20, 21  and  22  into  eq  13; 

AQ  =  C,  (1-CJ  S  A  I  (1  -w)  /  +  o  (7,^  Q  a  -  7;*  )  (1  -C,)  +  aCJT:^  -  ) 

+  bUp^L,  (e.  -  Co)  +  bUp^L,  7(7,  -  7^)  +  RC^  (7p  -  T„)~Kg  37/3z|^  .  (23) 

The  heat  transferred  by  the  different  processes  changes  the  heat  content  of  the  snowpack.  A  net 
energy  influx  will  reduce  the  energy  deficit  of  the  pack  and  if  the  pack  energy  deficit  is  erased  and 
the  pack  becomes  isothermal  at  0°C,  any  additional  energy  influx  will  cause  melting.  On  the  other 
hand  a  net  outflow  of  heat  will  cause  freezing  of  the  free  water  in  the  pack  and/or  an  increase  in 
the  energy  deficit.  These  phenomena  may  be  expressed  by 

AQ  =  xp,  C,  A7,  +  Lfp^  AW  (24) 

where  x  is  the  snowpack  depth,  AT  is  the  change  in  snow  temperature  and  AW  is  the  change  in 
water  in  the  snowpack. 

The  water  in  the  snowpack  comes  either  from  melting  or  rain.  This  water  may  freeze  or  remain 
as  liquid  within  the  pack.  When  the  pack  becomes  isothermal  at  0°C  and  becomes  saturated,  any 
additional  water  from  melt  or  rain  becomes  outflow  at  the  bottom  of  the  pack. 

Substituting  eq  24  into  eq  23  gives  the  heat  balance  in  the  snowpack: 

xp,  C, A7,  +  LfP„ AW  =  C,  (1  -C^) S-A  L{l-(J)I  +  aiT^C^a-T*)(\-C^) 

+  -  K)  +  bUp„L^  (c,  -  Co)  +  bUp^L^  y(T,  -  7^)  +  RC^  (7p  -  7o) 

-KgbTIbz]^.  (25) 
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The  water  equivalent  of  the  snowpack  should  be  held  constant  during  the  calculation  of  the  energy 
balance.  When  a  mass  change  takes  place  due  to  either  precipitation  or  meltwater  outflow,  the 
water  equivalent  of  the  pack  must  be  adjusted  before  calculating  the  heat  balance;  only  the  heat 
content  of  the  mass  change  is  included  in  the  calculation. 

In  most  situations  some  of  the  terms  in  eq  25  may  be  neglected.  The  heat  transfer  at  the  air/snow 
interface  is  usually  much  larger  than  that  at  the  soil/snow  interface  if  the  ground  is  frozen,  so  the 
latter  term  may  be  neglected.  When  the  ground  is  unfrozen  the  heat  flux  is  usually  small  enough 
to  be  neglected,  except  for  shallow  snowpacks  when  short-wave  radiation  penetrates  to  the  ground 
and  causes  heating  at  the  interface.  The  topography,  vegetation  and  time  during  the  snow  season 
will  determine  which  of  the  other  terms  become  ne^gible. 

Heat  transfer  equation 

The  heat  transfer  within  the  snowpack  must  be  calculated  to  determine  the  change  in  heat  stor¬ 
age.  The  heat  transfer  mechanisms  within  the  pack  are  conduction,  radiant  energy  transfer,  and 
phase  change. 

The  mathematical  expressions  for  these  processes  are  related  to  the  physical  characteristics  of 
the  snow,  such  as  density,  grain  size  and  degree  of  intergranular  bonding.  Because  the  latter  two 
quantities  are  hard  to  measure,  density  is  used  to  characterize  the  snowpack. 

Since  heat  is  gained  or  lost  by  a  snowpack  through  the  upper  and  lower  boundaries,  the  temper¬ 
ature  gradients  within  the  pack  are  vertical.  The  conduction  process  is  therefore  a  one-dimensional 
process  and  may  be  expressed  by  the  Fourier  conduction  equation: 

/:,32  7;/3z2=CiP.37’,/3r  (26) 

where  is  the  effective  thermal  conductivity  and  t  is  time. 

The  effective  conductivity  is  used  to  account  for  the  combined  effects  of  heat  transfer  by  con¬ 
duction  and  phase  changes.  Rewriting  this  equation: 

d^TJdz^=K^dTJdt  (27) 

where  =  C-^pJK^  . 

Schwerdtfeger  (1963)  proposed  a  relationship  between  the  diffusivity  and  density: 

K^=2K.J{3p,-p;)C^  (28) 

where  K^  is  the  thermal  conductivity  of  ice  and  Pj  and  p,  are  the  densities  of  ice  and  snow. 

In  a  natural  snowpack  the  density  increases  with  increasing  distance  from  the  upper  boundary. 
Incorporating  this  depth  dependence  into  eq  27  and  substituting  eq  28  into  the  revised  conduction 
equation  yields 

32 7’,/3z2  =  dTJdt  +  ^3--^  3p,/3z  dTJdz  .  (29) 

If  the  heat  transfer  by  the  latent  heat  processes  in  the  snow  is  accounted  for  by  the  use  of  the 
effective  thermal  conductivity,  this  equation  can  be  used  to  develop  a  temperature  profile  of  the 
snowpack. 

Snowpack  density 

Snow  is  a  dynamic  medium.  For  mobility  purposes  it  is  necessary  to  know  the  density  of  the 
snowpack  at  any  time  in  its  history.  The  metamorphic  process  is  one  of  densification,  with  four 
processes  involved:  compaction,  equitemperature  or  destructive  metamorphism,  temperature 
gradient  metamorphism,  and  melt  metamorphism. 
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Tlic  coinpaction  of  snow  layers  in  seasonal  snowpacks  has  been  observed  by  several  investigatois 
(U.S.  Army  Corps  of  Kngineers  l‘)56,  Kojiina  1967).  Kojinia  derived  an  expression  relating  the 
change  in  density  of  a  layer  to  the  weight  of  the  overlying  layers.  This  expression,  combined  with 
a  relationship  among  snow  temperature,  viscosity  and  density  derived  by  .Mellor  (1964)  yields 

1/p,  ^pJ^t=WC^  exp  [-0.08(27.^16-7)1  exp  j-Cj  pj  (30) 

where 

W  =  weight  of  snow  compressing  the  layer  of  interest,  expressed  as  water  equivalent 
C^  =  empirically  determined  constant  for  the  fractional  increase  in  density  at  0°C  and 
Ps  =  0 

Cj  =  empirically  detennined  constant. 

The  constants  depend  on  the  snow  type.  Equation  30  may  be  used  for  all  conditions  except  low- 
density  new  snow,  depth  hoar  or  windpacked  snow. 

Destructive  metamorphism  is  a  process  in  which  the  water  molecules  rearrange  themselves  on 
the  snow  crystals  so  that  the  free  surface  energy  is  reduced.  It  is  by  this  process  that  the  various 
snowflakes,  with  their  relatively  complicated  shapes  and  large  surface-to-volume  ratios,  are  meta¬ 
morphosed  into  rounded,  smooth  ice  grains.  Tliis  is  how  a  newly  fallen  snow  layer  is  incorporated 
into  the  existing  snowpack,  with  a  marked  settling  of  the  layer  and  an  increase  in  density.  Anderson 
(1976)  has  developed  an  expression  for  densification  by  this  process; 


1 /p,  3p,/ar  =  C3  exp  [-C4  (273.16- T,)]  exp  [^6(p,-pj)]  (31) 

where 

C3  =  fractional  settling  rate  at  0°C  for  densities  tess  than  Pj  (hr-') 

C4  =  empirically  determined  constant 

46  =  constant  (cm^/gr)  necessary  to  reduce  the  settling  rate  when  Pj  <  p, 

Pj  =  empirically  determined  density. 


This  equation  incorporates  the  increasing  rate  of  metamorphism  with  increasing  temperature  and 
the  preference  for  destructive  metamorphism  in  the  new  snow  layers,  which  tend  to  be  less  dense. 

Constructive  metamorphism  takes  place  when  there  is  a  temperature  gradient  in  the  snowpack. 
Water  vapor  moves  from  the  warmer  parts  of  the  snow  to  the  cooler  parts,  producing  a  density 
gradient.  The  increase  in  temperature  also  makes  the  pore  air  unsaturated.  Some  of  the  vapor  may 
not  be  transferred  but  may  be  used  to  keep  the  pore  air  saturated.  If  the  effective  diffusion  coef¬ 
ficient  73g  for  water  vapor  is  a  function  of  the  local  temperature  in  the  pack,  the  rate  of  change  of 
density  with  temperature  is  (Anderson  1976) 


*  dr 


(32) 


where  Cis  the  water  vapor  concentration. 

Melt  metamorphism  is  important  during  the  spring  when  the  pack  becomes  isothermal,  during 
rains,  and  when  the  air  or  ground  temperatures  become  greater  than  O’C.  In  a  low-density  fresh 
snow  layer,  water  increases  destructive  metamorphism  and  the  settling  rate.  Freeze-thaw  cycles 
in  general  cause  compaction  and  thus  a  density  increase  in  the  affected  layer.  Anderson  (1976) 
introduced  a  multiplicative  parameter  into  the  density  change  equation  for  destructive  metamor¬ 
phism.  This  factor  is  greater  than  one  and  has  to  be  determined  empirically  whenever  water  is 
present. 
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Implementation  of  the  model 

Calculation  of  the  energy  balance  (eq  25)  and  the  conduction  (eq  29)  requires  data  describing 
both  the  terrain  and  the  meteorology.  A  calculation  is  performed  to  determine  the  snow  conditions 
at  a  particular  location.  These  snow  conditions  are  then  assumed  to  be  the  same  over  the  surround¬ 
ing  area  or  terrain  unit.  The  extent  of  a  terrain  unit  is  determined  by  the  homogeneity  of  the  slope 
aspect,  elevation  and  vegetation  of  the  area.  The  meteorological  data  needed  are  cloud  cover, 
temperature,  wind  speed,  dewpoint  or  relative  humidity,  precipitation,  and  snow  and  ground  tem¬ 
peratures.  The  following  section  is  an  explanation  of  the  terms  in  the  energy  balance  and  conduc¬ 
tion  equations  with  comments  on  the  constants  and  the  method  of  calculation  where  further  expla¬ 
nation  is  needed. 

Tire  short-wave  radiation  reaching  the  snow  surface  consists  of  two  components,  direct  and 
diffuse.  The  direct  radiation  is  that  fraction  of  the  solar  radiation  transmitted  througli  tlie  atmos¬ 
phere  without  being  scattered  or  absorbed.  The  diffuse  radiation  reaches  the  snow  surface  after 
being  scattered  in  the  atmosphere  and/or  reflected  by  the  surrounding  terrain.  Diffuse  radiation 
reaches  the  snow  surface  from  all  directions  but  is  most  intense  in  the  direction  of  the  sun.  Al¬ 
though  diffuse  radiation  may  be  a  significant  fraction  of  the  total  short-wave  energy  flux,  it  is 
not  easy  to  express  mathematically  and  is  often  ignored  in  energy  balance  calculations.  In  this 
model  only  the  direct  beam  is  calculated  (in  the  first  term  on  the  right  side  of  eq  25),  althougli 
some  of  the  forward  scattered  energy  may  be  accounted  for  by  increasing  the  cloud  cover  factor 

The  amount  of  solar  radiation  transmitted  through  the  atmosphere  to  the  snow  surface  depends 
on  the  latitude,  time  of  year,  time  of  day,  slope  and  orientation  of  the  surface,  atmospheric 
transparency,  cloud  cover  and  vegetative  cover.  The  seasonal,  latitude  and  surface  orientation 
effects  can  be  calculated  using  a  method  formulated  by  Gamier  and  Ohmura  (1968).  If  the  mer¬ 
idian  of  solar  noon  is  used  as  a  reference,  then  the  motion  of  the  sun  is  north  and  south  along  this 
meridian,  and  the  unit  vector  S,  representing  the  sun’s  position,  may  be  specified  by  the  declination 
5  (Fig.  1): 

5  =(0,  cos5,sin5) .  (33) 


Zenith  g 


Figure  1.  Geometry  us^  to  calculate  the  insolation 
on  a  swface  (shaded).  X  =  unit  normal  to  the  sur¬ 
face,  S  =  position  vector  of  the  sun,  6  =  declination, 
A  =  azimuth  and  Z  =  zenith  angle. 
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The  orientation  of  the  snow  surface  is  described  with  respect  to  an  altitude-azimuth  coordinate 
system,  with  the  axes  in  the  local  horizon  plane  and  going  through  the  zenith.  The  unit  normal  to 
the  surface  is  expressed  in  terms  of  the  slope  an^e  and  azimuth: 

X  =  (-coSi4  sinZ,  siny4  sinZ,  cosZ)  (34) 

where  A  is  the  azimuth  angle  (north  =  0°,  east  =  90°,  etc.)  and  Z  is  the  zenith  angle. 

This  vector  is  then  transformed  to  the  global  coordinate  system: 

X"  =  {(-sin^sin//  cosy4  sinZ  +  cos//siny4  sinZ 
+  cos  0  sin  //  cos  Z),  (-sin  0  cos  H  cos  A  sin  Z 
+  sin  //  sin  y4  sin  Z  +  cos  0  cos  H  cos  Z),  (cos  0  cos  A  sin  Z 

+  sin0cosZ)]  (35) 

where  0  is  the  latitude  and  H  is  the  hour  angle  measured  from  noon  (//  is  negative  in  the  forenoon 
and  positive  in  the  afternoon). 

The  reduction  in  insolation  due  to  the  angle  at  which  the  radiation  strikes  the  slope  is  equal  to 
cosine  of  the  angle  between  the  sun’s  direction  and  the  line  normal  to  the  slope: 

^  -S  =  cosH 

=  (-sin  0  cos//  cos /I  sinZ  +  sin //sin /I  sinZ 
+  cos  0  cos  H  cosZ)  cos  6 

+  (cos0cosy4  sinZ  +  sin0cosZ)sin6,  (36) 

In  addition  to  the  effects  of  geometry,  the  energy  reaching  the  snowpack  is  reduced  by  atmos¬ 
pheric  absorption  and  scattering.  Garnier  and  Ohmura  (1968)  presented  a  simple  relationship  due 
to  Haltner  and  Martin  (1957),  by  which  the  energy  /incident  on  a  horizontal  surface  may  be  cal¬ 
culated  from  the  atmospheric  transmissivity  and  the  solar  constant: 

/=/oP'"  (37) 


where 

Iq  =  solar  constant  (1.35x10^  J/m^  s) 
P-  mean  zenith  path  transmissivity 
m  =  optical  air  mass. 


The  optical  air  mass  is  the  ratio  of  the  atmospheric  path  length  traveled  by  the  sun’s  rays  to  the 
atmospheric  depth  at  the  zenith.  This  may  be  obtained  from  tables  (list  1963,  Table  137).  It 
may  also  be  calculated  using  the  solar  elevation  angle  a  (Eagleson  1970): 


sino  =sin6  sin0  +  cos5  cos0  =  cos 


and  an  approximation  given  by  Choudhury  (1979a): 


(38) 


m  = 


1 

sina 


for  a  >  30° 


and 
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m  =  -A-  -  2.8/a2  for  a  <  30°  . 
sin  a 


(39) 


For  a  clear  sky  the  factor  in  the  first  term  representing  incident  short-wave  radiation  in  the 
energy  balance  equation  can  now  be  written  as 


S-A-L-I 


^0  J//, 


P'”  cos  HdH 


(40) 


where  the  limits  of  integration  are  the  sun’s  hour  angle  at  the  beginning  and  end  of  the  energy  bal¬ 
ance  calculation  interval.  The  limiting  values  of  and  H2  are  the  local  sunrise  and  sunset,  respec¬ 
tively.  For  flat  terrain  the  local  sunrise  and  sunset  may  be  obtained  from  The  American  Ephemeris 
and  Nautical  Almanac*  or  calculated  by  solving  eq  38  with  a  =  0.  For  areas  where  the  topography 
changes  the  local  sunrise  and  sunset  times,  the  times  found  for  flat  terrain  may  be  substituted  into 


sinZj  =  -cos6  sin///cosa  (41) 

which  is  solved  for  Z^,  the  azimuth  of  the  sun.  With  the  azimuths  of  sunrise  and  sunset,  the  alti¬ 
tudes  of  the  topographical  features  at  these  azimuths  is  then  substituted  into  eq  38  to  calculate  the 
local  sunrise  or  sunset  times. 
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Optical  Air  Mosi 

Figure  2.  Generalized  transmission  through 
clouds  as  a  function  of  cloud  type  and  air 
mass.  (From  List  1963.) 


The  effects  of  cloud  cover  on  the  radiation 
budget  are  expressed  as  averages  for  each  in¬ 
terval  of  integration: 

C^  =  \~N  +  NT  (42) 

where  is  the  fractional  cloud  cover  in  tenths 
and  r  is  a  cloud  transmission  coefficient. 

The  cloud  transmission  coefficient  de¬ 
pends  on  the  cloud  type,  thickness  and 
hei^t;  no  general  relationship  exists  be¬ 
tween  these  quantities.  Figure  2  shows 
transmission  as  a  function  of  cloud  type  and 
optical  air  mass  (List  1963).  These  measure¬ 
ments  do  not  separate  direct  and  diffuse  ra¬ 
diation,  partially  compensating  for  the  lack 
of  a  separate  diffuse  radiation  term. 

The  albedo  w  depends  on  many  factors, 
including  size,  shap)e  and  structural  arrange¬ 
ments  of  the  snow  crystals,  snowpack  water 
content,  snowpack  age,  sun  angle,  cloud 
cover,  and  amount  of  foreign  material  on 
the  surface.  Newly  fallen  snow  generally 
has  a  high  albedo,  0.8  or  higher.  This  de¬ 
creases  as  the  snow  ages  and  can  attain  val¬ 
ues  as  low  as  0.4  for  pure  snow.  The  albedo 
decrease  is  correlated  with  grain  size  increase, 
densification,  temperature  and  time.  Temp¬ 
erature  is  the  independent  variable  that  af¬ 
fects  albedo  the  most,  because  the  rate  of 
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snowpock) 


Age  of  Snow  (doysi 

a.  The  albedo  as  a  function  of  the  number  of  days  since  the  last  snowfall. 
The  prairie  albedo  shows  effects  of  the  underlying  surface  (Male  and  Gray 
1975).  The  mountain  snowpack  curves  are  from  U.S.  Army  Corps  of  Engi 
neers  (1956). 


Summotion  of  Temperotures  I  F) 

b.  The  albedo  as  a  function  of  the  sum  of  the  maximum  daily  temperatures 
since  die  last  snowfall  (U.S.  Army  Corps  of  Engineers  1956). 


Figure  3.  Decay  of  snowpack  albedo. 


Nat  Radiotion  to  Snow  (W/cm^) 


snow  metamorphism  depends  on  it.  Figure  3  illustrates  the  decrease  in  albedo  as  functions  of 
temperature  and  time.  The  difference  in  albedo  between  the  accumulation  and  ablation  seasons 
illustrates  how  the  higher  temperatures  during  the  late  winter  and  spring  increase  the  rate  of  snow 
metamorphism  from  that  during  the  accumulation  season.  The  curve  for  the  shallow  prairie  snow- 
pack  illustrates  the  effects  of  the  underlying  surface  on  the  albedo  as  bare  patches  appear  and  the 
remaining  snow  becomes  thin  enough  for  the  underlying  surface  to  affect  the  albedo. 

Indexing  the  albedo  to  temperature  provides  a  general  guide  for  clean,  deep  snowcover  (Fig.  3b). 
Where  the  snow  is  contaminated  with  other  material  (soot  or  dust,  for  example),  the  initial  albedo 
is  lower  than  that  of  pure  snow.  As  the  contaminated  snow  melts,  more  low-refl>  ctivity  material 
is  exposed,  so  that  the  albedo  decays  faster  during  the  ablation  season  than  during  the  accumula¬ 
tion  season.  The  relationship  for  the  albedo  of  clean  snow  as  a  function  of  temperature  is 


u=  140.15  -  14.4  In  (7)  (43) 

where  Tis  the  sum  of  the  daily  maximum  temperatures  since  the  last  snowfall. 

The  two  long-wave  radiation  exchange  terms  are  a  (1-C„)  [aC^  ^  ^  ]  for  the  exchange  be¬ 
tween  the  snowpack  and  the  sky  and  aC^[T^  -  TJ]  for  the  excnange  between  Uie  snowpack  and 
the  vegetation  canopy.  The  vegetative  cover  factor  ranges  from  0  for  no  canopy  to  1  for  a  com¬ 
plete  cover,  such  as  a  dense  evergreen  canopy. 

Figure  4  illustrates  the  effect  of  the  canopy  on  the  radiation  budget  of  the  snowpack  on  a  clear 
day.  The  canopy  has  the  greatest  effect  on  the  incoming  short-wave  radiation.  As  the  canopy 
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Figure  4  Effect  of  canopy  density  on  radiant 
energy  exchanges  at  the  snow  surface  (Jeffrey 
1968). 


density  increases,  the  loss  in  incoming  short¬ 
wave  energy  is  partially  compensated  for  by 
the  decrease  in  long-wave  energy  loss  to  the 
sky  and  the  increase  in  long-wave  energy  re¬ 
ceived  from  the  canopy.  The  net  radiational 
energy  received  by  the  snowpack  is  the  alge¬ 
braic  sum  of  three  components.  The  net 
energy  received  by  the  snowpack  has  a  min¬ 
imum  at  about  20%  canopy  density  and  then 
increases  with  increasing  canopy  density. 

The  amount  of  net  radiation  received  under 
a  canopy  of  density  1  is  about  half  that  for 
the  same  conditions  without  a  canopy. 

Figure  4  represents  the  radiation  balance 
during  the  day  when  the  canopy  is  relatively 
warm.  The  canopy  temperature  varies  di- 
urnally,  being  influenced  most  by  the  air 
temperature  and  insolation.  The  insolation 
may  raise  the  temperature  of  the  unshaded 
parts  of  the  canopy  several  degrees  above 
the  ambient  air  temperature,  but  the  effect 
of  this  on  the  radiation  to  the  snow  is  small, 
so  the  canopy  temperature  may  be  approx¬ 
imated  by  the  air  temperature. 

The  temperatures  and  temperature  gra¬ 
dients  within  the  snowpack  vary  on  a  daily 
basis  due  to  the  diurnal  cycle  and  on  a  longer 
time  scale  associated  with  weather  patterns 
The  incident  solar  radiation  warms  the  snow 
indirectly  through  higher  daytime  air  tem¬ 
peratures  and  directly  through  radiation 
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absorption  in  the  top  50  cm  of  the  snowpack.  The  diurnal  variation  of  air  temperature  is  also  ex¬ 
hibited  by  the  snowpack.  The  amplitude  of  the  variations  in  the  snowpack  decreases  with  increas¬ 
ing  depth  and  the  phase  lag  incrertes  with  depth,  with  most  of  the  effects  damped  out  at  a  depth  of 
15-20  cm  (Quick  1967,  LaChapelle  and  Armstrong  1977).  The  variations  in  internal  snowpack 
temperatures  due  to  weather  patterns  also  exhibit  damping  and  phase  lag  effects  but  are  irregular. 

Although  the  snowpack  is  relatively  transparent  to  short-wave  radiation,  it  is  highly  absorbing 
in  the  infrared,  so  that  the  temperature  of  the  top  centimeter  governs  the  infrared  radiation  char¬ 
acteristics.  The  temperature  of  this  surface  layer  can  be  considered  to  be  the  same  as  the  air  tem¬ 
pera'  ure  for  calculating  the  long-wave  radiation  balance  except  when  the  air  temperature  is  above 
0°C.  For  air  temperatures  greater  than  this  the  snow  temperature  is  set  at  0°C. 

The  terms  for  the  latent  and  sensible  heat  exchange,  bUp^L^  -  e^)  and  bUp^L^y  -  Tg), 
depend  on  the  conditions  in  the  turbulent  boundary  layer,  the  first  two  or  three  meters  above  the 
snow  surface.  The  turbulence  in  this  layer  is  determined  by  the  wind,  the  insolation,  and  local  con¬ 
ditions  such  as  surface  roughness,  topography  and  vegetation.  These  terms  are  evaluated  by  meas¬ 
uring  the  air  temperature  and  humidity  gradients  and  the  surface  wind  at  the  point  of  interest. 

DEVELOPING  INPUT  DATA 

The  model  presented  here  has  been  developed  from  components  that  require  extensive  input 
data.  The  feasibility  of  using  this  model  for  predicting  the  snowpack  characteristics  that  affect 
mobility  depends  on  whether  the  input  data  can  be  extrapolated  from  measurements  made  at  loca¬ 
tions  other  than  where  the  model  is  to  be  applied. 

Several  techniques  are  used  for  extrapolating  measurements  from  one  location  to  another. 
Anderson  (1973)  developed  a  technique  for  extrapolating  the  air  temperature  at  a  point  from 
measurements  at  several  locations  bracketing  that  point.  Sli^tly  different  methods  are  adapted 
to  mountainous  areas,  where  elevation  changes  have  a  significant  effect  on  the  temperature,  and 
to  low  relief  areas,  where  elevation  does  not  affect  the  temperature  significantly.  These  methods, 
possibly  refined  by  using  estimates  of  insolation  and  surface  conditions,  can  be  used  to  estimate 
the  air  and  snow  temperatures.  The  snow  albedo  may  also  be  estimated  using  Figure  3b. 

Kuz'min  (1961)  developed  a  relationship  for  estimating  the  humidity  as  a  function  of  altitude. 
This  may  apply  to  mountainous  regions,  but  the  influence  of  local  environment  and  surface  condi¬ 
tions  is  also  very  important;  a  general  relationship  describing  these  influences  has  yet  to  be  de¬ 
veloped. 

Topographic  and  reconnaissance  maps  are  readily  available,  so  that  the  insolation  and  exposure 
can  be  found  easily.  The  recoimaissance  maps  will  probably  contain  the  necessary  information  on 
the  vegetative  canopy.  A  quantitative  estimate  of  the  canopy  is  therefore  available,  but  its  relation¬ 
ship  to  the  snowpack  is  difficult  to  quantify. 

Meiman  (1968)  tabulated  the  results  of  24  investigations  of  forest  cover  and  accumulation. 

These  investigations  covered  a  wide  variety  of  conditions  and  had  widely  varying  results.  Some  of 
the  studies  using  multiple-factor  regression  analysis  showed  that  accumulation  is  inversely  corre¬ 
lated  to  canopy  density.  This  agrees  with  the  results  shown  in  Figure  5.  Canopy  density  nullifies 
the  effects  of  wind  (Fig.  6)  and  small-scale  topography,  producing  a  more  uniform  snowfall  distri¬ 
bution  under  the  canopy. 

Forest  cover  affects  the  ablation  phase  of  the  snowcover  by  reducing  the  net  radiation  at  the 
snow  surface,  reducing  the  melt  rate.  Hendrie  and  Price  (1979)  found  that  the  turbulent  energy 
exchange  for  a  deciduous  forest  was  reduced  to  the  point  where  it  can  be  neglected,  and  although 
the  net  solar  radiation  is  reduced,  it  can  be  used  to  calculate  the  energy  balance.  FfoUiott  et  al. 
(1965)  studied  the  effects  of  ponderosa  pine  stands  on  the  accumulation  and  ablation  of  snow  in 
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Figure  5.  Effect  of  vegetative  canopy  on  the  accumulation  of  snow  under  the 
canopy  (U.S.  Army  Corps  of  Engineers  1956). 


Wind  (mph) 

Figure  6.  Effect  of  vegetation  on  the  wind  profile  for  two  tree  heights.  (After  Jeffrey 
1 968  and  Izumi  1968.) 


adjacent  clearings.  Figure  7  illustrates  the  effect  of  the  stand  on  the  distribution  of  snowcover  as 
a  function  of  distance  from  tlie  forest  edge.  For  the  prevailing  winds  the  clearing  was  in  the  lee  of 
the  forested  area  so  that  the  accumulation  pattern  was  due  to  wind  transport  of  snow  from  the 
forest  to  the  forest  edge. 

During  the  ablation  phase  the  melt-rate  distribution  in  Figure  7  was  governed  by  radiative  ex¬ 
change.  Tlie  forest  was  to  the  southwest  of  the  clearing,  reducing  the  short-wave  radiation  inci¬ 
dent  on  the  snowpack.  This  produces  a  situation  analogous  to  that  in  Figure  4  (to  the  left  of  the 
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Oistonce  from  Canopy  (units  of  canopy  height) 

Figure  7.  Modification  of  mowpack  distribution  near  the 
edge  of  a  stand  of  trees.  The  top  graph  shows  the  typ¬ 
ical  snowpack  distribution  during  the  accumulation  phase, 
while  the  bottom  graph  shows  the  distribution  in  melt 
rate  due  to  the  canopy.  (After  Ffolliott  et  al.  1965.) 

minimum  in  the  all-wave  gain  curve),  where  increasing  distance  from  the  canopy  corresponds  to 
decreasing  canopy  density.  These  figures  illustrate  how  the  forest  cover  and  wind  combine  to 
affect  the  snowpack  accumulation  in  Figure  7,  and  forest  cover  and  the  aspect  of  the  forest  edge 
affect  the  solar  radiation  and  the  resulting  melt. 

Vegetation  on  prairies  and  steppes  affects  snow  accumulation  in  an  entirely  different  manner. 

In  these  regions,  where  the  ground  is  relatively  flat,  the  vegetation  (generally  either  grass  or  crop 
stubble)  affects  the  surface  roughness.  The  variation  in  surface  roughness,  combined  with  the  pre¬ 
valent  winds  on  the  prairies,  produces  a  wide  variation  in  snow  depths.  Moderate  winds  and  a 
relatively  uniform  vegetation  cover  produces  a  snowpack  with  a  depth  equal  to  the  vegetation 
height,  because  additional  snow  is  blown  away.  Hi^  winds  can  strip  away  the  accumulation  or 
prevent  the  pack  from  forming.  Variations  in  vegetation  height  can  produce  an  even  snowcover 
or  a  cover  of  drifts  and  bare  ground. 

The  effects  of  the  wind  are  probably  the  most  difficult  to  quantify.  As  shown  in  the  above 
examples,  the  effects  of  each  of  the  environmental  factors  such  as  vegetation  are  modified  by  wind. 
The  wind  also  plays  an  important  role  in  the  heat  exchange  processes  at  the  air/snow  interface. 

The  transfer  of  sensible  and  latent  heat  are  governed  by  turbulent  exchange  in  the  boundary  layer. 
This  turbulence  is  controlled  in  part  by  the  wind  speed  and  wind  speed  profile.  Male  and  Granger 
(1979)  showed  the  relationships  among  the  stability  parameter  and  the  gradients  of  wind  and  tem¬ 
perature,  from  which  the  turbulent  energy  terms  may  be  derived.  Scott  (1964)  presented  a  relation¬ 
ship  between  the  large-scale  “gradient  wind”  and  surface  roughness,  from  which  the  wind  profile 
in  the  region  near  the  surface  layer  may  be  obtained.  If  these  relationships  can  be  refined  through 
further  research,  it  may  be  possible  to  predict  the  turbulent  energy  exchange  terms. 

The  wind  also  has  an  effect  on  the  snowcover  density.  In  the  absence  of  wind  the  density  of 
newfallen  snow  depends  on  the  crystal  structure  of  the  snowflakes  and  varies  between  10  and  190 
kg/m^  (Diamond  and  Lowry  1953).  The  lighter  snows  are  formed  of  plate  and  dendritic  snowflakes. 
When  these  flakes  are  fractured  by  wind-induced  saltation,  they  can  pack  into  a  denser  snowcover. 
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with  densities  approaching  210  kg/m^ ,  depending  upon  the  wind  speed.  Bilello  (1969)  showed  the 
combined  effects  of  wind  and  temperature  on  density  to  be 

p  =  0.152  -  0.0031  7"  + 0.019  V  (44) 

where 

p  =  average  seasonal  snowcover  density 
T  =  average  seasonal  air  temperature 
V  =  average  seasonal  wind  speed. 

Although  this  relationship  applies  only  to  seasonal  averages,  it  illustrates  the  importance  of  wind 
speed.  There  have  been  many  observations  ofsnowpack  metamorphosis  with  time  and  tempera¬ 
ture,  but  there  do  not  appear  to  be  any  general  relationships  that  can  be  used  to  predict  the  model 
input  parameters.  The  components  of  the  model  that  provide  for  the  increase  in  density  due  to 
metamorphosis  and  compaction  require  a  number  of  empirical  constants.  Methods  of  predicting 
these  constants  and  other  relationships  of  metamorphism  have  yet  to  be  developed. 

In  addition  to  these  limitations  the  model  is  limited  because  the  various  model  components 
have  been  developed  for  relatively  deep  snowpacks.  For  snowpacks  less  than  about  20  cm  the 
model  results  don’t  apply  ;  this  would  not  affect  the  predictions  for  most  vehicles,  because  snow 
that  shallow  would  be  completely  compacted  under  the  weight  of  the  vehicle.  Then  the  snow 
acts  as  a  low  adhesion  layer,  and  its  other  properties  are  negligible. 

CONCLUSIONS 

The  modeling  elements  put  forth  in  this  study  can  provide  a  basis  for  a  model  to  predict  snow- 
pack  parameters  that  effect  mobility.  In  general  the  extensive  snowpack  characteristics  can  be 
predicted  with  more  accuracy  and  from  fewer  measurements  than  the  intensive  characteristics. 

Two  gross  snow  parameters,  depth  and  density,  are  probably  the  most  significant  for  mobility 
use.  Predicting  the  depth  requires  the  development  of  a  suitable  accumulation  model.  This  model 
could  be  based  on  the  weather  forecast,  with  elements  to  adjust  the  accumulations  for  specific 
terrain  and  vegetation  features.  Once  the  pack  depth  is  determined,  the  density  can  be  determined 
if  a  method  can  be  developed  to  remotely  determine  the  empirical  model  parameters. 

The  model  elements  that  have  been  presented  are  a  starting  point  from  which  to  build  the 
model.  To  develop  the  model  further  will  require  effort  in  two  areas:  1)  refinement  and  defin¬ 
ition  of  terrain  elements  over  which  the  snowpack  can  be  considered  uniform,  and  2)  development 
of  either  simple  techniques  to  measure  parameters  such  as  the  vapor  pressure  and  snow  densification 
constants,  or  models  to  eliminate  the  need  to  measure  these  parameters. 

Ultimately  the  feasibility  of  the  model  depends  on  the  accuracy  with  which  the  desired  param¬ 
eters  can  be  predicted.  Since  the  model  inputs  are  a  weather  forecast,  topographical  data  and 
parameters  extrapolated  from  one  or  more  sets  of  remote  measurements,  the  uncertainty  of  the 
inputs  can  be  large.  The  utility  of  the  model  will  depend  on  the  reduction  of  these  uncertainties 
to  an  acceptable  level. 

LITERATURE  CITED 

Adams,  W.P.  (1976)  Areal  differentiation  of  snow  cover  in  east  central  Ontario.  Water  Resources 
Research,  12(6):  1226-1234. 


Adams,  W.P.  and  D.R.  Barr  (1976)  The  evolution  of  snowcover  in  east  central  Ontario  with  special 
reference  to  the  International  Field  Year  for  the  Great  Lakes  (IFYGL)and  to  stratigraphy.  Pro¬ 
ceedings  of  the  33rd  Annual  Eastern  Snow  Conference,  pp.  1-21 . 

Anderson,  E.A.  (1973)  National  Weather  Service  River  Forecast  Svstem-Snow  accumulation  and 
ablation  model.  NOAA  Technical  Memorandum  NWS  HYDRO-17,  National  Weather  Service, 

Slver  Spring,  Maryland. 

Anderson,  E.A.  (1976)  A  point  energy  and  mass  balance  model  of  a  snow  cover.  NOAA  Technical 
Report  NWS  19,  Office  of  Hydrology,  National  Weather  Service,  Silver  Spring,  Maryland. 

Anderson,  H.W.  and  AJ.  West  (1965)  Snow  accumulation  and  melt  in  relation  to  terrain  in  wet 
and  dry  years.  Proceedings  of  the  33rd  Annual  Western  Snow  Conference,  pp.  73-S2. 

Baker,  M.B.,  Jr.  and  D.R.  Carder  (1977)  Comparative  evaluations  of  four  snow  modeb.  Proceed¬ 
ings  of  the  45th  Annual  Western  Snow  Conference,  pp.  58-62. 

Banon,  H.  (1978)  A  Markov  model  of  snow  depth.  Massachusetts  Institute  of  Technology,  Depart¬ 
ment  of  Civil  Engineering,  Research  Report  R78-1. 

Klello,  MA.  (1957)  A  survey  of  arctic  snow-cover  properties  as  related  to  climatic  conditions. 
SIPRE  Research  Report  39.  AD-158  191. 

Klello,  M.A.  (1969)  Relationships  between  climate  and  regional  variations  in  snow-cover  density 
in  North  America.  CRREL  Research  Report  267.  AD-700  990. 

BileOo,  M.A.  (1971)  Frozen  precipitation:  Its  frequency  and  associated  temperatures.  Proceed¬ 
ings  of  the  28th  Annual  Eastern  Snow  Conference,  pp.  68-80. 

Qioudliiuy,  BJ.  (1979a)  Thermal  conditions  in  a  snowpack.  I.  The  energy  balance  compmients. 
Report  prepared  by  Computer  Sciences  Corp.  for  National  Aeronautics  and  Space  Administration 
under  Contract  NAS5-24350. 

Cboudhury,  BJ,  (1979b)  Thermal  conditions  in  a  snowpack.  11.  Temperature  dbtribution  due  to 
tolar  heating.  Report  prepared  by  Computer  Sciences  Corp.  for  NASA  under  contract  NAS5- 
24350. 

Otoudhury,  BJ,  (1979c)  Snow  albedo  variation  with  solar  altitude  and  fractional  cloud  cover. 
Report  prepared  by  Computer  Sciences  Corp.  for  NASA  under  contract  NAS5-24350. 

Cuiry,  G.E.  and  A.S.  Mann  (1965)  Estimating  precipitation  on  a  remote  headwater  area  of  western 
Alberta.  Proceedings  of  the  33rd  Annual  Western  Snow  Conference,  pp.  58-66. 

Diamond,  M.  and  W.P.  Lowiy  (1 953)  Correlation  of  the  density  of  new  snow  with  700  mb  temper¬ 
ature.  SIPRE  Research  Paper  1 . 

Fames,  Pi.  (1971)  Mountain  precipitation  and  hydrology  for  snow  surveys.  Proceedings  of  the 

39th  Annual  Western  Snow  Conference,  pp.  44-49. 

t*^e»oa,t,S.  (1970)  Dynamic  Hydrology.  New  York:  McGraw-Hill. 

FfoUiott,  P.F.,  E.A.  Hanson  and  AJ).  Zander  (1965)  Snow  in  natural  openings  and  adjacent 
ponderosa  pine  stands  on  the  Beaver  Creek  watersheds.  U.S.  Forest  Service  Research  Note  RM-53. 
Flebiger,  G.  (1976)  The  influence  of  relief  and  vegetation  on  snow  deposition  in  the  Alps  with 
special  regard  to  the  effects  on  skiing  areas.  Proceedings  of  the  33rd  Annual  Eastern  Snow  Con¬ 
ference,  pp.  156-162, 

Galway,  J.G.  (1952)  The  effect  of  altitude  and  latitude  on  New  England  snowfall  estimates.  Pro¬ 
ceedings  of  the  9th  Annual  Eastern  Snow  Conference,  pp.  57-61, 

Gamier,  BJ.  and  A.  Ohmura  (1968)  A  method  of  calculating  the  direct  shortwave  radiation  income 
of  slopes.  Journal  of  Applied  Meteorology,  7(5):  796-800. 

Gary,  H.L.  (1975)  Airflow  patterns  and  snow  accumulation  in  a  forest  clearing.  Proceedings  of 
the  43rd  Annual  Western  Snow  Conference,  pp,  106-113. 

Gray,  D.M.  (1979)  &iow  accumulation  and  dbtribution.  Proceedings  of  a  Meeting  on  Modeling 
of  Snow  Cover  Runoff,  26-28  September  1978,  Hanover,  New  Hampshire.  CRREL,  pp.  3-33. 
Haltner,  G  J.  and  R.F.L.  Martin  (1957)  Dynamical  and  Physical  Meteorology.  New  York:  McGraw- 
Hill. 


Hendrick,  R.L.  and  R  J.  DeAngelis  (1976)  Seasonal  snow  accumulation,  melt  and  water  input— A 
New  England  model.  Journal  of  Applied  Meteorology,  15(7):  717-727. 

Hendrie,  L.K.  and  A.G.  Price  (1979)  Energy  balance  and  snow  melt  in  a  deciduous  forest.  Pro¬ 
ceedings  of  a  Meeting  on  Modeling  of  Snow  Cover  Runoff,  26-28  September  1978,  Hanover,  New 
Hampshire.  CRREL,  pp.  2 1 1  -22 1 . 

Isyumov,  N.  and  A.G.  Davenport  (1974)  A  probabilistic  approach  to  the  prediction  of  snow  loads. 
Canadian  Journal  of  Gvil  Engineering,  1 :  28-49. 

Isyumov,  N.  and  M.  Mikitiuk  (1976)  Climatology  of  snowfall  and  related  meteorological  variables 
with  application  to  roof  load  specifications.  Proceedings  of  the  33rd  Annual  Eastern  Snow  Con¬ 
ference,  pp.  41-69. 

bumi,  Y.  (1968)  Kansas  1968  field  data  program.  Air  Force  Cambridge  Laboratories,  Hanscom 
Air  Force  Base,  Bedford,  Massachusetts. 

Jeffrey,  W.W.  (1968)  Snow  hydrology  in  the  forest  environment.  Snow  Hydrology,  Proceedings 
of  the  Workshop  Seminar.  Canadian  National  Committee,  International  Hydrological  Decade,  pp. 
1-19. 

Kind,  R.J.  (1976)  A  critical  examination  of  the  requirements  for  model  simulation  of  wind-induced 
erosion/deposition  phenomena  such  as  snow  drifting.  Atmospheric  Environment,  10(3);  219-227. 
Kojima,  K.  (1967)  Densification  of  seasonal  snow  cover.  Physics  of  Snow  and  Ice-Proceedings  of 
the  International  Conference  on  Low  Temperature  Science-Sapporo,  Vol.  1 ,  Part  2.  The  Instit- 
tute  of  Low  Temperature  Science,  Hokkaido  University,  Sapporo,  Japan,  pp.  929-952. 

Kuz'min,  P.P.  (1961)  Protsess  tayaniya  skezhnogo  prokrova.  Gidrometeorolgicheskoe  Izdatef stvo, 
Leningrad  (Melting  of  Snow  Cover.  Israel  Program  for  Scientific  Translations,  TT  71-50095, 1972). 
LaChapelle,  E.R.  and  R.L.  Armstrong  (1977)  Temperature  patterns  in  an  alpine  snow  cover  and 
their  influence  on  snow  metamorphism.  Institute  of  Artie  and  Alpine  Research,  University  of 
Colorado,  Boulder,  Colorado. 

Leaf,  C.F.  and  G.E.  Brink  (1973)  Computer  simulation  of  snowmelt  within  a  Colorado  subalpine 
watershed,  USDA  Forest  Service  Research  Paper  R1 1-99,  Rocky  Mountain  Forest  and  Range  Ex¬ 
periment  Station,  U.S.  Forest  Service,  Ft.  Collins,  Colorado. 

List,  RJ.  (1963)  Smithsonian  Meteorological  Tables.  Smithsonian  Institution,  Washington,  D.C. 
Male,  D.H.  and  D.M.  Gray  (1975)  Problems  in  developing  a  physically  based  snowmelt  model. 
Canadian  Journal  of  Gvil  Engineering,  2:  474-488. 

Male,  D.H.  and  R.J.  Granger  (1979)  Energy  and  mass  fluxes  at  the  snow  surface  in  a  prairie  en¬ 
vironment.  Proceedings  of  a  Meeting  on  Modeling  of  Snow  Cover  Runoff,  26-28  September  1978, 
Hanover,  New  Hampshire.  CRREL,  pp.  101-124. 

McKay,  G.A.  (1968)  Problems  of  measuring  and  evolution  of  snowcover.  In  Snow  Hydrology, 
Proceedings  of  the  Workshop  Seminar.  Canadian  National  Committee,  International  Hydrological 
Decade,  pp.  49-62. 

McKay,  G.A.  and  B.F.  Findlay  (1971)  Variation  of  snow  resources  with  climate  and  vegetation 
in  Canada.  Proceedings  of  the  39th  Annual  Western  Snow  Conference,  pp.  17-25. 

McKay,  D.C.  and  G.W.  Thurtell  (1978)  Measurement  of  the  energy  fluxes  involved  in  the  energy 
budget  of  a  snow  cover.  Journal  of  Applied  Meteorology,  17(3):  339-349. 

McHirtland,  J.T.,  A.B.  Super  and  V.L.  Mitdiell  (1971)  Snowpack  accumulation  in  relation  to  terrain 
and  meteorological  factors  in  southwest  Montana.  Report  prepared  for  Army  Research  Office— 
Durham,  Durham,  North  Carolina,  by  Montana  State  University. 

Meiman,  J.R.  (1968)  Snow  accumulation  related  to  elevation,  aspect  and  forest  canopy.  Snow 
Hydrology,  Proceedings  of  the  Workshop  Seminar.  Canadian  National  Committee,  International 
Hydrological  Decade,  CRREL,  pp.  35-47. 

MeIlor,M.  (1964)  Properties  of  snow.  CRREL  Monograph  III-Al.  AD-611  023. 

Myers,  V.A.  (1962)  Airflow  on  the  windward  side  of  a  large  ridge.  Journal  of  Geophysical  Research, 
67(1):  4267-4291. 


Niloiy*,  U.  (1954)  Snow  OysJlaii:  Natural  and  Artificial.  Cambridge,  Massachusetts:  Harvard 
University  Press. 

hGeben,  B.C.  (1966)  A  technique  for  forecasting  the  rate  of  snowfall  at  Alta,  Utah.  Alta  Avalanche 
Study  Center,  Project  D,  Progress  Report  #1,  U.S.  Dept,  of  Agriculture,  Forest  Service,  Wasatch 
National  Forest. 

Obled,  C.  and  H.  Harder  (1979)  A  review  of  snow  melt  in  the  mountain  environment.  Proceedings 
of  a  Meeting  on  Modeling  of  Snow  Cover  Runoff ,  26-28  September,  1978,  Hanover,  New  Hamp¬ 
shire.  CRREL,pp.  179-204. 

Quick,  M.C.  (1967)  A  comparison  of  measured  and  theoretical  snowpack  temperatures.  Journal 
of  Hydrology,  5:  1-20. 

Rhea,  J.O.  and  LO.  Grant  (1974)  Topographic  influences  on  snowfall  patterns  in  mountainous 
tenain.  In  Advanced  Concepts  and  Techniques  in  the  Study  of  Snow  and  Ice  Resources.  Wash¬ 
ington,  D.C.:  National  Academy  of  Sciences,  pp.  182-192. 

Rhea,  J.O.  (1978a)  Orographic  precipitation  model  for  hydrometeorological  use.  Department  of 
Atmosi^eric  Sciences,  Colorado  State  University,  Atmospheric  Science  Paper  No.  287. 

Rhea,  J.O.  (1978b)  Performance  of  a  simple  orographic  precipitation  model  for  hydrological  ano 
climatological  studies.  Presented  at  ASCE  Spring  Convention,  Pittsburgh,  Pennsylvania. 
Reifsnyder,  W.E.  (1955) Wind  profiles  in  a  small  isolated  forest  stand.  Forest  Science,  1:  289-297 
Sage,  J.D.  (1976)  An  operational  model  for  hourly  snowfall.  Proceedings  of  the  33rd  Annual 
Eastern  Snow  Conference,  pp.  98-128. 

Scott,  R.F.  (1964)  Heat  exchange  at  the  ground  surface.  CRREL  Cold  Regions  Science  and  Engi¬ 
neering  Monograph  IIAl .  AD-449  934. 

Schwefdtfeger,  P.  (1963)  Theoretical  derivation  of  the  thermal  conductivity  and  diffusivity  of 
snow.  International  Association  of  Scientific  Hydrology,  Publication  No.  61,  pp.  75-81. 

Spreen,  W.C.  (1947)  A  determination  of  the  effect  of  topography  upon  precipitation.  Transac¬ 
tions  of  the  American  Geophysical  Union,  28(2);  285-290. 

Steppuhn,  H.  and  G.E.  Dyck  (1974)  Estimating  true  basin  snowcover.  In  Advanced  Concepts  and 
Techniques  in  the  Study  of  Snow  and  Ice  Resources.  Washington,  D.C.:  National  Academy  of 
Sciences,  pp.  314-328. 

Stoir,  D.  and  D.L  Golding  (1974)  A  preliminary  water  balance  evaluation  of  an  intensive  snow 
survey  in  a  mountainous  watershed.  In  Advanced  Concepts  and  Techniques  in  the  Study  of  Snow 
and  Ice  Resources.  Washington,  D.C.:  National  Academy  of  Sciences,  pp.  294-303. 

Tabler,  R.D.  (1975)  Predicting  profiles  of  snowdrifts  in  topographic  catchments.  Proceedings  of 
the  43rd  Annual  Western  Snow  Conference,  pp.  87-97. 

U.S.  Army  Corps  of  Engineen  (1 956)  Snow  Hydrology:  Summary  Report  of  Snow  Investigations. 
North  Pacific  Division,  Portland,  Oregon. 

Vance,  H.M.  and  B.L.  Waley  (1971)  &iow  frequency  analysis  for  Oregon  and  Utah.  Proceedings 
of  the  39th  Annual  Western  Snow  Conference,  pp.  34-38. 

Willen,  D.W.,  CJV.  Shumway  and  J.E.  Reid  (1971)  Simulation  of  daily  snow  water  equivalent  and 
melt.  Proceedings  of  the  39th  Annual  Western  Snow  Conference,  pp.  1-8. 

Willis,  N.O.  (1979)  Snow  on  the  Great  Plains.  Proceedings  of  a  Meeting  on  Modeling  of  Snow 
Cbver  Runoff,  26-28  September  1978,  Hanover,  New  Hampshte.  CRREL,  pp.  56-62. 


26 


A  facsimile  catalog  card  In  Library  of  Congress  MARC 
format  Is  reproduced  below. 


Berger,  Roger  H. 

Developing  a  model  for  predicting  snowpack  param¬ 
eters  affecting  vehicle  mobility  /  by  Roger  H.  Berger. 
Hanover,  N.H.:  U.S.  Cold  Regions  Research  and  Engi¬ 
neering  Laboratory;  Springfield,  Va.:  available  from 
National  Technical  Information  Service,  1983. 

V,  34  p.,  lllus.;  28  cm.  (  CRREL  Report  83-16.  ) 
Bibliography:  p.  23. 

1.  Mobility.  2.  Snow.  3.  Heat  balance. 

4.  Vehicle  mobility.  I.  United  States.  Army.  Corps 
of  Engineers.  II.  Cold  Regions  Research  and  Engi¬ 
neering  Laboratory,  Hanover,  N.H.  III.  Series: 

CRREL  Report  83-16. 


ruicsi  cuvcr  iujccis  aoidutjii  pnasc  ui  uxc  Mxuwyuvcr  oy  icuucjji^  uic  nci  lauuxuvii  ai  uic 

snow  surface,  reducing  the  melt  rate.  Hendrie  and  Price  (1979)  found  that  the  turbulent  energy 
exchange  for  a  deciduous  forest  was  reduced  to  the  point  where  it  can  be  neglected,  and  although 
the  net  solar  radiation  is  reduced,  it  can  be  used  to  calculate  the  energy  balance.  FfoUiott  et  al. 
(196S)  studied  the  effects  of  ponderosa  pine  stands  on  the  accumulation  and  ablation  of  snow  in 


